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ABSTRACT 

It is widely believed that the cosmological redshift is not a Doppler shift. However, Bunn & Hogg 
have recently pointed out that to settle properly this problem, one has to transport parallely the ve- 
locity four- vector of a distant galaxy to the observer's position. Performing such a transport along 
the null geodesic of photons arriving from the galaxy, they found that the cosmological redshift 
is purely kinematic. Here we argue that one should rather transport the velocity four-vector along 
the geodesic connecting the points of intersection of the world-lines of the galaxy and the ob- 
server with the hypersurface of constant cosmic time. We find that the resulting relation between 
the transported velocity and the redshift of arriving photons is not given by a relativistic Doppler 
formula. Instead, for small redshifts it coincides with the well known non-relativistic decompo- 
sition of the redshift into a Doppler (kinematic) component and a gravitational one. We perform 
such a decomposition for arbitrary large redshifts and derive a formula for the kinematic compo- 
nent of the cosmological redshift, valid for any FLRW cosmology. In particular, in a universe with 
f2 m = 0.24 and fi A = 0.76, a quasar at a redshift 6, at the time of emission of photons reaching 
us today had the recession velocity v — 0.997c. This can be contrasted with v = 0.96c, had the 
redshift been entirely kinematic. Thus, for recession velocities of such high-redshift sources, the 
effect of deceleration of the early Universe clearly prevails over the effect of its relatively recent 
acceleration. Last but not least, we show that the so-called proper recession velocities of galaxies, 
commonly used in cosmology, are in fact radial components of the galaxies' four-velocity vectors. 
As such, they can indeed attain superluminal values, but should not be regarded as real velocities. 

Key words: methods: analytical - cosmology: theory 
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1 INTRODUCTION 

A standard interpretation of the cosmological redshift in the frame- 
work of the Friedman-LemaTtre-Robertson- Walker (FLRW) models 
is that it is an effect of the expansion of the Universe. This inter- 
pretation is obviously correct since 1 + z = a(t )/a(t c ), where z 
is the value of the redshift, a(t) is the scale factor of the Universe 
and t e and t are respectively the times of emission and observa- 
tion of a sent photon. In semi-popular literature (e.g. Kaufmann & 
Freedman 1999; Franknoi, Morrison & Wolff 2004; Seeds 2007), but 
also in professional (e.g. Harrison 2000; Abramowicz et al. 2007), 
one can often find statements that distant galaxies are 'really' at rest 
and the observed redshift is caused by the 'expansion of space'. Ac- 
cording to other authors (e.g. Peacock 1999; Whiting 2004; Chodor- 
owski 2007a,b; Bunn & Hogg 2009, hereafter BH9), such statements 
are misleading and cause misunderstandings about the cosmological 
expansion. A presentation of the ongoing debate in the literature on 
this issue is beyond the scope of the present work; a (possibly non- 
exhaustive) list of papers includes Davis & Lineweaver (2001), Davis 
& Lineweaver (2004), Barnes et al. (2006), Francis et al. (2007), 
Lewis et al. (2007), Lewis et al. (2008), Gr0n & Elgar0y (2007), Pea- 
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cock (2008), Abramowicz et al. (2009), Chodorowski (2008), Cook 
& Burns (2009). 

On the other hand, there is broad agreement that the cosmolog- 
ical redshift is not a pure Doppler shift; the gravitational field must 
also generate a gravitational shift. Gravity can be neglected only lo- 
cally, in the local inertial frame (LIF) of an observer. It turns out that 
for small redshifts, the cosmological redshift can be decomposed into 
a Doppler shift and a Newtonian gravitational one (Bondi 1947). The 
latter is a shift induced by the Newtonian gravitational potential. Can 
the cosmological redshift be decomposed into a Doppler shift and a 
gravitational shift (not necessarily Newtonian) for an arbitrary value 
of the redshift? This is the question which we want to deal with in this 
Paper. Formally, the answer is no. There is no invariant definition of 
the recession velocity of a distant galaxy in general relativity (GR). 
This velocity is a relative velocity of the galaxy and the observer, 
and in curved spacetime there is no unique way to compare vectors 
at widely separated points. A natural way to define the recession ve- 
locity is to transport parallely the velocity four- vector of the distant 
galaxy to the observer, but the result will depend on the chosen path. 
(This is just the definition of curvature.) In practice, however, as a 
'preferred' path one can choose a geodesic connecting the galaxy and 
the observer. Moreover, in FLRW models there is a natural foliation 
of spacetime, into space-like hypersurfaces of constant cosmic time. 
In our Paper, as the geodesic we will adopt the geodesic lying on such 
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a hypersurface, i.e. connecting the points of intersection of the world- 
lines of the galaxy and the observer with the hypersurface of constant 
cosmic time. 

In a seminal study of the cosmological redshift, BH9, following 
Synge (1960) and Narlikar (1994), adopted another geodesic for the 
parallel transport: the null geodesic along which the photon is trav- 
elling from the source to the observer. This approach results in one 
'effective' velocity, while we think it is important to make a distinc- 
tion between the velocity at the time of emission and the velocity at 
the time of observation. These two velocities are obtained by trans- 
porting parallely the velocity four-vector of the source respectively 
on the hypersurface of constant t e and t Q . Not surprisingly, our result 
differs from that obtained by BH9. However, unlike theirs, ours cor- 
rectly reproduces the small-redshift decomposition of the cosmologi- 
cal redshift into a Doppler component and a gravitational component, 
mentioned above. Specifically, our decomposition coincides with that 
of Bondi (1947) for isotropic and homogeneous matter distribution. 

This Paper is organized as follows. In Section [2] we transport 
parallely the velocity four- vector of a distant galaxy to the observer. 
From the transported vector we calculate the recession velocity of 
the galaxy, which turns out to depend on the galaxy's comoving 
distance and the assumed background cosmological model. In Sec- 
tion [3] we find specific relations between the cosmological redshift 
and its Dopplerian component for two particularly simple cosmologi- 
cal models: the empty model and the Einstein-de Sitter model. In Sec- 
tion [4] using only the Principle of Equivalence, we derive the small- 
redshift decomposition of the cosmological redshift and find it to be 
identical with that obtained in Section [3] using generally relativistic 
approach. In Section [5] we calculate the recession velocities in the 
two models mentioned above, as well as for the currently favoured, 
flat non-zero A model with Q m = 0.24. We find also that the ve- 
locities are subluminal in all FLRW cosmological models and for all 
values of the redshift. A comparison of the results obtained in this 
Paper with some earlier works on the subject is given in Section [6] 
Summary is presented in Section|7] 
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Figure 1. Parallel transport of the velocity four- vector of a distant galaxy, U, 
from the galaxy to the observer along the geodesic on the hypersurface of 
constant cosmic time. The transport is plotted for two values of time: the ob- 
servation time, t (dashed horizontal line), and the emission time, t a (solid 
horizontal line), both in units of t a . The observer is in the origin of the coor- 
dinate system, x = 0, and the comoving coordinate of the emitting galaxy is 
(vertical line) in units of cH^ 1 . The transport is performed in comoving 
coordinates; in these coordinates all spatial components of the galaxy's veloc- 
ity four-vector are initially zero. During the transport towards the observer, the 
radial component of the four-velocity gradually increases. For better visibility 
and reasons explained in Subsection l5.2l in the plot this component is multi- 
plied by a(t), Slanted long-dashed lines show the direction of null vectors at 
the events (t , 0) and (t e , 0) in local inertial coordinates of the observer. 



2 PARALLEL TRANSPORT 

The equation for parallel transport of a vector U along a curve with a 
parameter A and a tangent vector U = {dx^/dX} is 



dU° dx» _ 

~dx +v ^l^ u -°' 



(1) 



Here, F"^ are the Christoffel symbols, related to the underlying met- 
ric by the equation 



rot 1 a~t f . \ 

m/3 — 2~9 \9tP,n + 9iti,P - 9nPn) ■ 



(2) 



The metric considered here is the Friedman-Lemaitre- 
Robertson-Walker (FLRW) metric. In the RW coordinates, it is 



ds 2 = c 2 dt 2 — a 2 (t)[dx 2 + RqS 2 (x / Ro)dij) 2 
Here, 

dtp 2 = dd 2 + sin 2 6 dab 2 , 



(3) 



(4) 



and kRg is the present curvature of the universe. For a closed, flat, 
and open universe, respectively, k = +1,0,-1, and the function 
S(x) equals sin(a;), x, and sinh(a;). Time t is the cosmic time, that 
is the proper time of all Fundamental Observers partaking in the ho- 
mogeneous and isotropic expansion of the cosmic substratum (fluid). 
The scale factor a(t) relates fixed, or comoving, coordinates, x, to 



physical, or proper, coordinates, r: r = ax. The RW coordinates are 

{x^} = {ct,x,e,cp}. 

The curve considered here is a radial geodesic on the hypersur- 
face of constant cosmic time, hence dx^ — S^dx, where <5t is the 
Kronecker delta. This gives 



dU a 
dx 



= -r^ u 



(5) 



The vector U is here the four-velocity of a distant galaxy, parallely- 
transported to the central observer (i.e. to the origin of the coordi- 
nate system). By definition, U a = dx a /ds, where x a (t) denotes the 
galaxy world-line as a function of the cosmic time. In the RW (co- 
moving) coordinates every galaxy has fixed 8, <f> and x, so ds — cdt 
and U a = dx a /d(ct). For these reasons, the initial value for U a is 
5" } . Specifically, the initial conditions for Equations l[5) are 



U a (x*) = 5 a , 



(6) 



where x» is the comoving radial coordinate of the emitting galaxy. 

We calculate F"^ using equations {2}-l[3} and obtain the follow- 
ing set of two linear equations for U° and U 1 : 



dU° 
dx 

dU 1 
dx 



a 2 (t)H(t) 



H{t) 



U°. 



(7) 



(8) 
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We remind that since t is fixed, so are a and H . The corresponding 
equations for U 2 and U 3 are simple to solve; with conditions l[6j, 
the result is U 2 (x) — U 3 (x) = 0. Differentiating Equation 10 with 
respect to r and using l[8j to eliminate dll 1 /dx yields 



d 2 U l> 
dx 2 



a 2 H 2 u0 



The solution is 



U°{x) = de 



iHx j c 



+ C 2 e" 



:i H x J c 



(9) 



(10) 



where Ci and C2 are the integration constants. Imposing initial con- 
ditions and noticing that aH — awe obtain 



U°(x) = 



1 a(x — x*) /c 



1 a(x — x*)/c 



-d(x — x* ) / c 



1 — dfx — x±) /c 
+ 2^ e 



(ID 

(12) 



!7» - 2a . 
hence 

[/°(a; = 0) = cosh (das«/c) , (13) 

t/ 1 (a; = 0) = a -1 sinh (da;* /c) . (14) 

We see that (gapW^U 13 ) \ x=a = 1, as it should. 

In order to identify the transported velocity four-vector of the 
distant galaxy with its recession velocity, we now have to transform 
the vector from the RW coordinates to the coordinates of the Local In- 
ertial Frame (LIF) of the central observer. A general coordinate trans- 
formation 1° = x a (x a ) transforms the (unprimed) components of 
the metric, g a p, to 



g a 'p> 



dx a dx 13 



dx a ' dx? 



7 9af> ■ 



(15) 



The following transformation of the RW coordinates: t' = t, x' = 
a(t)x, 9' = 9,4> = 0, yields g a ipi which in the limit x' — I tend 
to rj a rpi, i.e. to the Minkowski metric []] Thus, the primed coordinates 
are indeed the coordinates of the observer's LIF. A transformation of 
vector components, U a ' = {dx a> / dx a )U a , yields here U ' = U° 
and U 1 ' = a(t)U l , hence 



U° {x = 0) = cosh(da;,/c) . 
U 1 {x = 0) = sinh (dx»/c) . 



(16) 
(17) 

are properly normalized: 



Again, the components U' 
(v a >f)'U a 'u ')\ x , =Q = l. 

In the LIF of the observer, the radial component of the parallely- 
transported four-velocity of the distant galaxy is non-zero. We in- 
terpret this effect as a non-zero recession velocity of the galaxy, 
v. Quantitatively, U° = 7 and U 1 = /?7, where ft = v/c and 
7 = (1 — /3 2 )~ 1//2 . This yields immediately 

7 = cosh (da;* /c) (18) 
/?7 = sinh(da;*/c) , (19) 

henc^l 

P = tanh {ax»/c) . (20) 
Finally, since sinh y + cosh y — exp y, we obtain 



1 To check this, instead of using eq. j\5\ it is much simpler to note that x = 
x' /a(t') and to take its full differential. 

2 The quantity ax t /c plays here essentially a role of the so called 'velocity- 
parameter' of the Lorentz transformation. See e.g. Rindler (1977) and Schutz 
(1985). 



exp 



1-/3 



1/2 



= 1 



(21) 



where z D is the Dopplerian component of the (total) cosmological 
redshift of the galaxy, z. In section[3]we will explain how to use this 
equation to relate the cosmological redshift of a distant galaxy to its 
Dopplerian component. Equation J2U is the central formula of the 
present paper. 



3 REDSHIFTS FOR SPECIFIC MODELS 

Let us remind that in Equation d21t . x* stands for the comoving ra- 
dial RW coordinate of a galaxy. For simplicity, from now on we will 
denote this coordinate by x. In fact, x is not only a coordinate but 
also the comoving radial proper distance to the galaxy and can be 
expressed in terms of the redshift of the galaxy. For a given cos- 
mological model, the relation between x and z is determined by 
the Friedman-Lemaitre (FL) equations. Using this relation and Equa- 
tion d21t enables one to deduce the relation between z and z D in this 
model. In the following, we will find the relation between z and z D 
for two specific, particularly simple cosmological models. 

Formula J21t is valid for any (fixed) value of cosmic time. Two 
instants of time are of interest: the time of observation, or 'today', t a , 
and the time of emission, t e . The corresponding values of the reces- 
sion velocity are respectively v a — v(t Q ) and v c — v(t e ). 



3.1 Time of observation 

At the time of observation a\ to = Ho (we use the normalization 
a = 1, so that the comoving radial distance equals to the proper 
distance today), hence 



exp 



H x(z) 



1 + ZD ■ 



(22) 



Here, Zu is a Doppler shift due to the present value of the recession 
velocity, v Q - 



3.1.1 Empty universe 

In the empty model (fi m = Qa ~ 0), the comoving radial distance is 
x(z) = cHg 1 ln(l + 2), hence e ln (1+z) = 1 + z D , or 



z = Zu ■ 



(23) 



In an empty universe the origin of the cosmological redshift is thus 
entirely Dopplerian. This result is expected, since spacetime of an 
empty universe is the Minkowski spacetime, where redshift is simply 
a Doppler shift. 

Some cosmologists still believe that the cosmological redshift is 
not a Doppler shift even in the case of an empty universe. The origin 
of this erroneous belief is a wrong interpretation of the FLRW met- 
ric for an empty universe in the RW coordinates. In these coordinates, 
space (3-D hypersurface of constant cosmic time) does have non-zero 
curvature, but what matters here is the full 4-D curvature, which is 
zero. Indeed, a simple coordinate transformation transforms the met- 
ric to the standard Minkowskian form. What else could spacetime 
of an empty universe be? There is no matter so there is no gravity, 
and a complete relativistic description is provided by special relativ- 
ity. Indeed, Chodorowski (2007a) proved that in an empty universe 
z — z-d, using only specially relativistic concepts. More specifically, 
he showed that 1 + jz d = a(t ) / a{t e ), where a(t) = t/t is the scale 
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factor of the empty model. Another derivation of the redshift in the 
empty model will be presented in Section[4] 

There is some remaining subtlety here that deserves further ex- 
planation. In the specially-relativistic formula for a Doppler shift, 
z D = + v/c)/(l — v/c), the recession velocity v is an inertial 
velocity of a distant galaxy, call it v distant - The velocity that appears 
in Equation ( 121b and resulting l !23t is a 'local' velocity, w local, trans- 
ported parallely from the galaxy to the central observer using the RW 
coordinates. Are these two velocities equal? Redshift is an observ- 
able, so the answer is 'they must be', but let us understand, why. 

In Minkowski, i.e. flat, spacetime, components of a transported 
four- vector in inertial coordinates are constant (see Eq.[T}, so U a loca i 
(i.e., transported) is equal to (7 a distant . Here, double primes de- 
note the global inertial coordinates of the central observer, in which 
v distant is measured. This equality is no longer true for non-inertial 
coordinates (like RW), but parallel transport does not depend on 
the coordinate system used, so the transported vector is the same. 
We have transformed the RW components of the transported four- 
velocity, £/ Q local , to the components in the coordinates of the LIF of 
the observer, !7 a local , and obtained «i OC al' Since the global inertial 
frame of the central observer is an extension of its LIF (possible only 
in an empty universe), U a local and U a local are the components of the 
same vector in the same coordinate frame, so they must be equal. We 
thus have !7 a ' distant = £/ Q 'i' oca i = (7 a ' local , hence v distant = «local. 

Summing up, our definition of the recession velocity of a distant 
cosmological object correctly implies that in the case of an empty 
universe the origin of the cosmological redshift is purely Dopplerian. 



3.2 Time of emission 

At the time of emission, a | t e = (a/a)\ c a c = H c /(1 + z), where H c 
is the Hubble constant at the time of emission. Using this equality in 
Formula J2U we obtain 



exp 



H{z)x(z) 



c(l + z) 



1 +ZD- 



(26) 



Writing the above we have used the fact that there is a unique cor- 
respondence between t e and z, so H c can be expressed in terms of 
the latter. A general formula for the Hubble constant as a function of 
redshift is 



H 2 (z) = Hi [fi A + Om(l + zf - (fi - 1)(1 + zf] 



(27) 



Here, D. = fi m + 17a- In Equation (26J, 2 D is a Doppler shift due to 
the recession velocity at the time of emission, v c . 



3.2.1 Empty universe 

From Equation i27l we obtain H(z) — Ho(l + z), hence 
H(z)x(z)/[c(l + z)] = In (1 + z). This yields again 



z u ■ 



(28) 



Equations l |28t and l |23t involve two different Doppler shifts: respec- 
tively at the time of emission, 2 D _ C , and at the time of observation, 
2 D o . These equations are consistent because in the empty model the 
expansion of the universe is kinematic and recession velocities are 
constant in time. This is not true for any other cosmological model. 



3.1.2 Einstein-de Sitter universe 

In the Einstein-de Sitter (EdS) model (fi m = 1 and f2 A = 0), the 
comoving radial distance is x 0) = 2cHo 1 [1-(1 + 2)" 1/2 ], hence 

2[l-(l + z)- 1/2 ] =ln(l + z D ). (24) 

The particle horizon for the EdS model is a; (2 = 00) = 2cff Q _1 . 
We see that the Dopplerian component of the infinite cosmological 
redshift of a source located at the horizon sphere is 2 D = e 2 — 1. 
This finite value corresponds to the recession velocity of the source 
not at the time of emission of photons (i.e., t e ~ 0), but today. In 
sub-subsection [3~2~2l we will see that at t c = 0, this velocity was c. 
The difference between the two velocities is due to deceleration of 
the EdS model. 



3.1.3 Small redshifts 

A second-order expansion in redshift of the radial distance valid for 
any cosmology (under the assumption that dark energy is in the form 
of the cosmological constant) is x(z) ~ cH^ 1 [2 — (1 + qo)z /2], 
where qo = fi m /2 — SI a- (It is easy to check that second-order ex- 
pansions of the exact formulas for the empty and EdS models are 
consistent with this general expansion.) From Equation d21t . after ex- 
pansion we obtain z D — 2 — go 2 2 /2, or 

z = 2 D + ^ 2 2 + O (4) • (25) 

In Section [4] we will see that the above result is identical to that ob- 
tained using the Principle of Equivalence (PoE). 



3.2.2 Einstein-de Sitter universe 

In the EdS model H(z) = H {1 + zf 2 , hence 

2 D = exp{2[(l + 2) 1/2 -l]} -1. (29) 

Infinite value of 2 corresponds to infinite value of 2 D , so at the time 
of the Big-Bang (corresponding to 2 = 00), the horizon was re- 
ceding from the observer with the velocity of light. Let us decom- 
pose the cosmological redshift into a Doppler shift and the 'rest', 
2 = 2 D + 'rest'. From Equation l |29t we see that for high redshifts, 
the Dopplerian component of the cosmological redshift is (much) big- 
ger than the total value of the latter. Therefore, in the EdS model the 
'rest' must in fact be negative, i.e., it is a blueshift. We will show later 
that for small redshifts, the 'rest' can be identified with a gravita- 
tional shift. It will be demonstrated that a gravitational shift is indeed 
a blueshift for all matter-only models. A homogeneous and isotropic 
distribution of ordinary matter creates a gravitational potential well 
with the observer at the center and the work performed by photons 
climbing down this well is negative. 



3.2.3 Small redshifts 

From Equation (27} we have H(z) = H [l + (1 + q )z + 0(z 2 )}, 
hence H(z) x(z)/[c(l + z)\ ~ 2 + (qo — l)2 2 /2. Inserted in Equa- 
tion (26}, this yields 

2 = 2 D - 22. z 2 + O (zl) . (30) 

Again, in Section [4] we will see that the above result is identical to 
that obtained using the PoE. 
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4 PRINCIPLE OF EQUIVALENCE 

In Section [5] we studied the cosmological redshift using generally- 
relativistic concepts. In this Section we will restrict our analysis to 
small values of the cosmological redshift; this will allow us to use the 
Principle of Equivalence (PoE). It is well known that the PoE enables 
one to study local effects of light propagation in a weak gravitational 
field. Specifically, one can replace a static reference frame in this field 
by a freely falling frame in empty space. From the PoE it follows that 
the latter is locally inertial, with the laws of light propagation be- 
ing given by special relativity. The classical examples of phenomena 
studied in this way are the gravitational redshift and the bending of 
light rays. In some instances it is possible to calculate correctly non- 
local effects by summing up infinitesimal contributions from adjacent 
local inertial frames, freely falling along light ray trajectory. For ex- 
ample, the (small) gravitational redshift of a photon sent from the 
surface of a non-relativistic star (like, e.g., the Sun) to infinity can be 
calculated in this way. In some other instances it is not possible; a 
classical example is the bending of light rays. The calculation based 
on the PoE gives only half of the value of bending predicted by GR 
in the weak-field approximation (but see Will 2006 and references 
therein). 

It is well known that the cosmological redshift can be inter- 
preted as 'an accumulation of the infinitesimal Doppler shifts caused 
by photons passing between fundamental observers separated by a 
small distance' (Peacock 1999; see also Peebles 1993, Padmanab- 
han 1993). The derivation is very simple but let us recall it here. 
Consider two neighbouring comoving observers, located on the tra- 
jectory of a light ray, separated by an infinitesimal proper distance 
81 = a(t)8x, where a(t) is the scale factor and 8x is their comoving 
separation. The velocity of the second observer relative to the first is 
Sv = (d/dt)Sl = aSx. The equation for the radial light propaga- 
tion in the RW metric is cSt = aSx, where St is the light travel time 
between the two observers, so Sv — cSa/a. This small recession ve- 
locity of the second observer causes a non-relativistic Doppler shift 
of the passing ray, <5A/A = 5v/c — 5a/ a. (Reference frames of the 
FOs are locally inertial.) Integrating the relation between A and a we 
obtain 1 + z = A /A e = a /a e . 

The above derivation resembles a derivation using the PoE, 
which will be presented below. There is, however, an important dif- 
ference. In the above we have applied the generally relativistic RW 
metric and the law of light propagation in it. The RW coordinates are 
local inertial coordinates of (freely falling) Fundamental Observers 
and the cosmic time t is their local proper time. This has enabled us 
to integrate the effect (along the ray's trajectory) trivially and exactly. 
Not surprisingly, then, the result is correct for an arbitrary value of 
the redshift. In an application of the PoE, the gravitational field of the 
system under consideration is described in terms of non-relativistic 
Newtonian dynamics, so the field must be weak. In the FLRW cosmo- 
logical models the gravitational field is weak only at small distances 
from an observer, more specifically much smaller than the Hubble ra- 
dius. Therefore, the PoE can give correct predictions only for small 
values of the cosmological redshift. Still, we think that calculations 
using the PoE are worth performing. The reasons are as follows. 

First, the aim of our Paper is to decompose the cosmological 
redshift into a global Doppler shift (i.e., due to the relative velocity 
of the emitter and the observer) and the 'rest'. In the above calcula- 
tion of the redshift we summed up contributions from infinitesimal 
relative velocities of neighbouring FOs, but this did not help us to ex- 
press the answer in terms of the relative velocity of the emitter and 
the observer. As stated already many times, for large cosmological 
distances the recession velocity a priori is not well defined. In this 



Paper we proposed its generally relativistic definition. This defini- 
tion is however quite formal, while for small redshifts, velocity is a 
well defined pre-relativistic concept. Therefore, the identification of 
a global Doppler component of the small cosmological redshift will 
be simple and unambiguous. Second, in previous sections, using GR, 
we have shown that the cosmological redshift, of any value, can be 
indeed decomposed into a Doppler shift and the 'rest'. In this section, 
using the PoE, we will similarly decompose the redshift into a global 
Doppler shift and the rest and will show that the 'rest' is the gravita- 
tional redshift in the Newtonian sense. In other words, we will prove 
that for small redshifts the rest is a classical Newtonian gravitational 
shift. 

Among all the FLRW models there is a toy model which can 
be described using special relativity only. This is the empty, or Milne 
model (Milne 1933). At first sight, the relation of this model to the 
PoE is problematic, since in an empty universe there is no gravity. 
However, the relevance of the model will be seen later. Spacetime 
of an empty universe is the Minkowski spacetime. Let us denote the 
global inertial coordinates of an observer by T, X, Y and Z. In the 
standard spherical coordinates, R, O, <E>, the metric is 

2 



ds z = dT A 



dR z 



R\d®* + sin ed-T 



(3D 



Like above, we will calculate the cosmological redshift summing 
up the infinitesimal Doppler shifts caused by photons passing be- 
tween neighbouring Fundamental Observers, located along the pho- 
tons' path. We have seen that such a summation can be easily done in 
the RW coordinates, for any FLRW model. However, it turns out that 
in the empty model, it is possible to perform this summation using 
also the Minkowskian coordinates T and R. 

In the Milne model the recession velocities of Fundamental Ob- 
servers are constant; for example, in the rest-frame O of the central 
observer v(T,R) — R/T — H(T)R — const. Let us consider a 
Fundamental Observer located at a distance R from the central ob- 
server at time T. In his rest-frame O', he observes a small Doppler 
shift of a photon sent (in fact, transmitted) by a nearby observer lo- 
cated at R + 5R. Due to the Hubble expansion, the latter observer has 
a small recession velocity Sv' = 5R' /T' — H(T')5R' relative to the 
former. Like in the RW coordinates (locally coinciding with the co- 
ordinates T' and R'), the relative wavelength shift is <5A/A = Sv' jc. 
However, in the empty model we can use the special-relativistic law 
of composition of velocities to relate the velocity Sv' to the veloc- 
ity Sv, measured in the observer's frame O. This law can be applied, 
because the relative velocity of the frames O and O', v(T, R), is con- 
stant. We have 



Sv' 



v(T,R + 5R) - v(T, R) 
l-v{T,R + 8R)v(T, R)/c 2 

hence we obtain 



Sv(T, R) 



[v{T,R)/c\ 



SX 

T 



8v(T, R)/c 



(32) 



(33) 



l-[v(T,R)/cY 

Let us denote the distance of the source at the time of emission by 
D e . Our present task is to integrate this equation along the photon's 
trajectory, c(T — T c ) = D c — R. (The photon travels towards the 
central observer.) However, it can be easily checked that since the 
velocity of any FO is constant, Sv(T,R) — H(T)8R equals to 
Sv e — H(T C )8R C . Therefore, the integration simplifies to an inte- 
gration over R c at fixed T — T c , 



dX 

T 



(H e /c) dR c 



(34) 



/ v - ./z,.i-(*W<0 2 ' 

(The minus sign accounts for the fact that for a photon travelling from 
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the source to the observer R decreases, so 8R is negative.) Defining 

/3 = H c R c /c and j3 e = H c D c /c, we obtain 



T 



5e 



d/3 
1 - p 2 



(35) 



The left-hand-side integral is equal to ln(l + z). The integration of 
the right-hand-side is elementary and yields ln[(l + / 8 G )/(1 — /3 e )]/2. 
Hence, 



1 + 2 



1/2 



1 + ZD- 



(36) 



In sum, we have accumulated infinitesimal Doppler shifts using the 
inertial coordinates of the central observer and expressed the total 
redshift in terms of the recession velocity of the source. The relation 
between the cosmological redshift and the recession velocity turned 
out to be given by the relativistic Doppler formula, in agreement with 
generally-relativistic Equation ( 128b . Thus, as already stated earlier, 
in an empty universe the cosmological redshift is entirely a global 
Doppler shift. A similar derivation of the redshift in the empty model 
using Minkowskian coordinates with an explicit integration along the 
photon's trajectory was presented in Appendix A and Section 2 of 
Chodorowski (2007a). 

Now we proceed to derive, using the PoE, the relation between 
the recession velocity of a nearby source and its redshift, for non- 
empty universes. Why was the exercise with the Milne model use- 
ful for this purpose? High redshifts correspond to relativistic veloc- 
ities and in the empty model their effects manifest as the 1 — f3 2 
factor in the integrand on the RHS of Equation d35 1 . Expansion of 
this integrand for small velocities and subsequent integration yields 
p c + (j8|). Therefore, in our non-relativistic calculations we will be 
allowed to retain terms up to second order in f3 c . (We will see later 
that the leading-order contribution to the velocity-redshift relation 
coming from Newtonian gravity is already of the second order in /3 C .) 
Like previously, we denote Newtonian coordinates of the central ob- 
server by T and R. They are now inertial only locally, because there 
is a gravitational field. We 'transform this field away' using local in- 
ertial frames of FOs. Contrary to the empty model, FOs do not move 
now with constant speeds. According to the Newtonian interpretation 
of the FL equations (Milne & Mc Crea 1934), the change of speed of 
any FO located on a sphere centered on another observer is induced 
by the gravitational force exerted on the FO by the mass encompassed 
within this sphere. (This is Newton's theorem for an isotropic and ho- 
mogeneous universe. The speed of the FO is, of course, measured rel- 
ative to the central observer). The effects of the cosmological constant 
A can be incorporated easily because A results in a force proportional 
to distance, similarly to the classical gravitational force for a homo- 
geneous mass distribution. The FO located at time T at a distance R 
from the central observer has velocity related to his velocity at the 
time of emission of a photon, u e , in the following way: 



v(T, R) = v c + g{T c , R C )AT + 0(AT 2 



(37) 



Here, T c is the time of emission, R c = R(T e ), AT = T - T c 
and g(T e , R c ) is the initial gravitational acceleration of the FO. Writ- 
ing the above formula we have approximated changing acceleration 
of the FO by its initial value. (The result of this approximation for 
v(T, R) is already of second order in AT.) It is straightforward to 
calculate g(T, R) from Newton's theorem and to express it in terms 
of the cosmological parameters. However, there is a faster way of do- 
ing this: by differentiating the Hubble law in Newtonian coordinates, 
v = H(T)R, with respect to time. The result is 

g(T,R) = -q(T)H 2 (T)R, (38) 



where q is the deceleration parameter. 

From Equation l|38ll, q(Tc.,R r ) = —q c H 2 R c ~ —qoH 2 R c . 
Inserting these expressions for v c and g in Equation l |37t , we obtain 
up to terms linear in AT 



v(T, R) = H C R C {1 - q H e AT) 



(39) 



(R ~ R c + H C R C AT). The velocity of a FO is described here 
in terms of T and his initial position, R c . A small relative veloc- 
ity of two neighbouring FOs at time T, 8v\ T , can be expressed as 
(dv/dR e )\ T 8R e = H c (l - q H c AT)5R c . This relative velocity 
causes a non-relativistic Doppler shift of the frequency of a passing 
photon. To obtain the total redshift we need to integrate these in- 
finitesimal Doppler shifts along the photon's path (i.e., null geodesic). 
Unlike in the empty model, here T and R are only locally inertial co- 
ordinates. However, their extent is quite limited, D c <C cH~ l , so 
H C AT equals to H C (D C ~ R)/c ~ H c (D a - R e )/c plus terms 
already of higher order in (D c — R c )/c. Therefore, d\/\ = 
ln(l + z) equals to 



1 



qo- 



H C {D C - R e 



dR c = /3 e 



1 «2 



(40) 



where, as before, /3 C 
1 



H c D c /c. Hence we obtain finally 



z — exp 



& - -<7o/3o 



^0t+O{pl). (41) 



Expansion of the Doppler formula d36>b to second order in /3 C gives 
z-d — Pc + jSc/2- Therefore, the first two terms on the RHS of 
Equation | |4U can be identified with a (global) Doppler shift due to 
the recession velocity of the source at the time of emission. Since 
the third term is of second order in /3 e , we can write qofi 2 /2 = 
goz D /2+0(z D ). In such a way we find out that Equation J4U is iden- 
tical to Equation d30t . the latter derived using a generally-relativistic 
approach. In Equation i25l , the kinematic component of the redshift 
is related by the Doppler formula to the recession velocity of the 
source at the time of observation. We recover this equation by noting 
(with some small labour) that Equation l |39t yields (3 C = /3 + qoPo> 
and using this equality in Equation J41t . 

The term ~qoP 2 /2 is a gravitational shift: the work, per unit- 
mass rest energy, performed by a photon travelling from the source 
to the observer. From Equation < I3 8b we have g ~ —qoH 2 R, so the 
gravitational potential ^(R) — qoH 2 R 2 /2 + const. We thus indeed 
obtain — $ e )/c 2 = —qoP 2 /2. Whether the observer is located in 
the centre of a potential well or on the top of a potential hill depends 
on the value of go- For all matter-only cosmological models go is 
positive, so the observer is in the centre of a potential well. Then a 
gravitational shift is negative, so it is in fact a blueshift. This reflects 
the fact that the photon climbs down a potential well so the performed 
work is negative. For an accelerating universe (like, as we currently 
believe, ours) go < and a gravitational shift is a redshift. 

Peacock (1999) derived Equation J4U in a much simpler way. 
He stressed that in the definition of the redshift, the wavelength at 
the time of emission is the wavelength measured in the source's 
rest-frame O' , A c '. In the observer's rest-frame O this wavelength 
is Doppler-shifted, A c /A c ' = 1 + z D . Next, on the path to the ob- 
server the photon undergoes a gravitational shift, A /A c = 1 + z G . 
Hence, 1 + z = A /A c ' = (1 + zd)(1 + z G ). This equality gives For- 
mula J41 1 because, as we have verified above, z G = —qof3 2 /2, and 
because z D z G is already of third order in /9 e . However, freely-falling 
observers employed in his calculation of z were not, strictly speak- 
ing, the freely-falling observers of the studied problem. In the RW 
models freely-falling observers are (by definition) FOs, partaking in 
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the Hubble expansion. Peacock employed the standard 'Einstein-lift' 
observers, starting to fall only when the photon passes the top of their 
cabin. Accounting for non-zero initial velocities of freely-falling ob- 
servers complicated our calculations, but enabled us to apply the PoE 
'from first principles'. We performed our calculations consistently in 
the same way as we did using global inertial coordinates in the empty 
model, and the RW coordinates in all FLRW models. 

Peacock's aim was to derive a formula for the angular diame- 
ter distance up to the terms quadratic in redshift. In order to obtain 
the correct result, he had to expand the relativistic Doppler formula 
up to the term quadratic in /3 C . Strictly speaking, this term is a rela- 
tivistic correction and in principle one is not allowed to account for 
it in the calculations using the PoE, where velocities are assumed to 
be non-relativistic. However, the justification for including this term 
can be found analyzing our calculations in the Milne model. We ob- 
tained ln(l + z u ) = /3 C + (D{$), or 1 + z D = exp[/3 c + O(^)], 
where O(pf) was the lowest-order relativistic correction. Therefore, 
solving for z a we were allowed to expand the exponential function 
up to a term quadratic in /3 e , still remaining in the non-relativistic 
regime. Since the exponential function has the same second-order ex- 
pansion as the Doppler formula, we correctly obtained the second- 
order Doppler relation, also in our non-relativistic calculations for 
non-empty models (see Eq.l4U. 

We stress that we did not assume that z — z u + z G . In our 
calculations both terms appeared naturally as a result of an applica- 
tion of the PoE. The first person to note that the cosmological red- 
shift, if small, can be decomposed into a Doppler shift and a grav- 
itational shift was Bondi (1947). He performed such a decomposi- 
tion in a spherically symmetric generalization of the RW models, the 
Lemaitre-Tolman-Bondi model. Our Formula J4U is a special case of 
his formula (52), in a sense that z G = —qofil/2 only for a homoge- 
neous distribution of cosmic matter and small distances. 

In the approach presented in this section, the cosmological red- 
shift of any value is generated by relative motions of FOs. For small 
values of the redshift the application of the PoE allowed us to de- 
compose the redshift into a Doppler component and a gravitational 
component|f] We have found that a Doppler component of the red- 
shift is obviously generated by a kinematic component of the motions, 
while a gravitational component of the redshift is generated by their 
varying-in-time component. Of course, this interpretation of a gravi- 
tational component is fully consistent with the standard interpretation 
of a gravitational shift, because this is a gravitational field that causes 
FOs to accelerate. 



5 RECESSION VELOCITIES 

Calculations in Section l(2j, based on our definition of recession ve- 
locity, led us to Equation 1201 . This equation is a general formula for 
the recession speed of a galaxy, valid for an arbitrary comoving po- 
sition of the galaxy and an arbitrary instant of time. The hyperbolic 
tangent on the RHS of this equation ensures that the recession veloc- 
ity of a source with an arbitrarily high redshift is subluminal (v < c). 
However, it is widely believed that recession velocities of sufficiently 
distant galaxies can be, and indeed are, superluminal (e.g. Davis & 
Lineweaver 2001; Davis & Lineweaver 2004; Gr0n & Elgar0y 2007; 
Lewis et al. 2007). In Subsection l l5.lt we will find a relation between 



our recession velocity and a commonly used estimator of the veloc- 
ity: the proper recession velocity. Basing on this relation we will show 
that the conception of the superluminal expansion of the Universe is 
wrong. In Subsection l !5.2t we will intercompare our recession veloc- 
ities of galaxies for some specific cosmological models. 



5.1 No superluminal velocities 

Let us remind that in the comoving coordinates {a; M } = {ct, x, 9, <f>}, 
the FLRW metric takes on the form given by Equation (3). For a par- 
ticle (or a galaxy) participating in the homogeneous and isotropic ex- 
pansion of the cosmic fluid, 9 = const, (f> = const and x = const. 
Hence, ds — cdt and 



^ = ^ = (1,0,0,0). 



(42) 



Let us change coordinates to {x^ } = {ct, r, 9, <f>\, where r — a(t)x. 
Spacetime interval ds is an invariant, so also in the primed coordi- 
nates ds — cdt. The four-velocity of the galaxy in these coordinates 



U a 



dx a 
cdt 



O.f.o,o) 



(43) 



where H (t) — a/ a is the Hubble constant. 

It is a standard practice in the literature to use as an estimator 
of the recession velocity of a galaxy its so-called proper recession 
velocity, 

dr 



dt 



(44) 



that is the derivative of the proper distance between the galaxy and 
the observer with respect to the proper cosmic time. Comparing equa- 
tion (IH with (Eg} we see that 



= cU x 



(45) 



or that Uprop/c is nothing more than the radial component of the 
galaxy's four-velocity in the proper coordinates {ct, r, 9, There- 
fore, the proper velocity can indeed attain superluminal valuesQ but 
should not be regarded as a recession velocity. As described earlier, 
to calculate the latter, one should 

A. Transport parallely the galaxy's velocity four-vector to the ob- 
server's position; 

B. Transform the components of the transported vector to the com- 
ponents in the locally inertial coordinates of the observer; 

C. Read off the value of relative velocity using the fact that in these 
coordinates, U 1 = /3~/(v). 

We stress that U a and U a given in Equations {42} and j43\ are just 
different components of the same four-vector U(t, r) in respectively 
comoving and proper coordinates. It has been noted earlier that in 
the vicinity of the observer, the proper coordinates reduce to her/his 
local inertial coordinates. Therefore, in this case the operation B. 
should not be necessary. Indeed, a straightforward (though somewhat 
more complicated) calculation of the parallel transport gives directly 
U (t, r = 0) = sinh (Hr/c), in agreement with Equation {T7). 

Note that for sufficiently distant objects, the transported (in 
proper coordinates) radial component of the four-velocity is even 
more superluminal than the initial one [sinh (Hr/c) compared to 



3 Another derivation of this decomposition was presented in Gr0n & Elgar0y 
(2007). 



4 For instance, in Fig.fTJthe present value of the proper velocity of the source 
is 1.24c. 
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Hr/c]E 

However, using the equation given in point C. of the above 
list, we get 



/3 OUT = tanh(v pr0 p/c) . 



(46) 



where w pr o P = Hr. Equation d46t . which coincides with equa- 
tion J20t , relates the proper recession velocity of a galaxy to the ve- 
locity which follows from our definition. Whatever the value of r 
(whence v plop ), v OUI remains smaller than c. 

The aim of this Paper is to derive a decomposition of the cosmo- 
logical redshift, so we have assumed that the distant galaxy lies on the 
past light cone of the observer. However, both Equations l[7]l-l[8j of 
the parallel transport and their initial conditions ((6) are completely in- 
dependent of this assumption. Consequently, in the solutions of these 
equations and in the resulting Equation J46b . the proper coordinate 
of a distant galaxy can be arbitrarily large; in particular, not limited 
by the radius of the particle horizon. As a corollary, not only the part 
of the cosmic substratum contained within the particle horizon, but 
entire cosmic substratum expands subluminally. Needless to say, the 
issue whether the Universe is open or closed does not change any- 
thing: even in an open universe, at any instant of time every particle 
has - though allowably arbitrarily large - finite proper coordinate. 

The only exception for subluminality of the expansion of the 
Universe is the moment of initial singularity. As then r — > 0, it 
is better to write Hr = ax. The early Universe is dominated by 
ultra-relativistic particles and radiation so its scale factor evolves as 
a(t) oc i 1//2 . Therefore, d oc i -1 / 2 , so it diverges at t = 0. The only 
point in space which in spite of this divergence is stationary at the 
Big-Bang is x = 0. This is obvious because x = is the origin of the 
coordinate system and all velocities are measured relative to it. How- 
ever, if the velocity of the central observer is measured by another 
observer, non-zero separation of the two observers and the fact that 
velocities are relative imply that the velocity of the central observer 
relative to all other observers is non-zero. (In order to prove it more 
formally, use the FLRW metric centered on any other point.) Thus, 
relative velocities of all FOs are non-zero and at the initial singular- 
ity all velocities attain the velocity of light. 

Summing up, the expansion of the Universe is never superlu- 
minal. A common misconception that the expansion is superluminal 
is based on the wrong identification of recession velocity with the 
'proper' recession velocity. 

This subject is interesting in its own right and worth of further 
study. Here we have presented it only briefly, to convince the reader 
that a decomposition of the cosmological redshift into a gravitational 
shift and a kinematic one for an arbitrary redshift is indeed possi- 
ble. This is so because for arbitrary redshifts, recession velocities of 
galaxies are subluminal. We plan to study in more detail the issue of 
(the lack of) superluminal expansion of the Universe on superhorizon 
scales in a follow-up work. 



5.2 Recession velocities for specific models 

Figure Q] shows parallel transport of the velocity four- vector of a dis- 
tant galaxy, U, from the galaxy to the observer along the geodesic 
on the hypersurface of constant cosmic time. The transport is plotted 
for two values of time: the observation time, t a (dashed horizontal 
line), and the emission time, t e (solid horizontal line), both in units 



5 In the model used to construct Fig.[T] the final value of the radial component 
transported on the hypersurface t = i is 1.58 compared to its initial value 
1.24. 



of t . The adopted cosmological model is a matter-only, open uni- 
verse with Sl m = 0.24. The observer is at the origin of the coordinate 
system, x — 0, and the comoving coordinate of the emitting galaxy 
is x, (vertical line) in units of cH^ 1 . The adopted redshift of the 
source at t e (determining the value of a;*) is equal to 3. The transport 
is performed in the comoving coordinates; in these coordinates all 
spatial components of the galaxy's velocity four-vector are initially 
zero. During the transport towards the observer, the radial component 
of the four- velocity gradually increases. In the plot this component 
is multiplied by a(t), to enable a physical interpretation of the trans- 
ported vector, 1/(0), in the framework of special relativity. Namely, 
at x — the scaled comoving components, (U°, all 1 ), are equal to 
the (non-scaled) components of U in the local inertial coordinates 
of the observer, (U® nl U^ n ). Slanted long-dashed lines show the di- 
rection of null vectors at the events (t o ,0) and (t e , 0) in the latter 
coordinates. We see in FigureQ]fhat E/£(0) [= aU 1 ^)] is greater for 
t e than for t . Thus, the recession velocity of the source at the time 
of emission was greater than it is at present. This conclusion is con- 
firmed by the directions of the vectors: both vectors are timelike (as 
required), but the direction of the vector at (t e , 0) is closer to the null 
direction than the direction of the corresponding vector at (t , 0) . It 
is expected, since a matter-only model constantly decelerates, so its 
expansion must have been faster in the past. 

In practice, it is more interesting to know the recession velocity 
of a source at the time it emitted the photons which are reaching the 
observer today. In other words, the velocities at the time of emission 
are more interesting than the velocities at the time of observation. 
(There is even no guarantee that the emitting galaxy exists today.) 

Combining Equation d20t with J2U and l !26t , we obtain 



/3 C = tanh 



H(z)x(z) 



c(l + z) 



(47) 



Again, Equation J47b assures that the recession velocity (at the time of 
emission) of a source with an arbitrarily high redshift is subluminal 
(v c < c). In the empty (Milne) model, from d47t we obtain /3 e = 
tanh[ln(l + z)\, or 



(1 + 2) 



(48) 



(l + 2) 2 + l ' 

in agreement with the Doppler formula. In the EdS model we obtain 

& = tanh{2[(l + 2) 1/2 -l]}. (49) 

In the currently favoured, flat model with non-zero A, Equation J27t 
yields H{z) = Hoyl — fi m + f2 m (l + z) 3 . The comoving dis- 
tance is x(z) = cH^JqII - f2 m + fim(l + y)' 3 ]~ 1/2 dy. 

Let us now consider a source located at a redshift, e.g., 6. For 
the Milne model, the corresponding recession velocity at the time of 
emission, v c M \ is 0.96c. For the EdS model, this velocity, vi , is 
about 0.997c. For the flat model with A, adopting Sl m = 0.24, the ve- 
locity, « e A) equals approximately to 0.991c. We interpret these values 
in the following way. A Milne universe is coasting, so « C M) = !) M) . 
An EdS universe decelerates, so it was expanding faster in the past. 
Therefore, it is logical that v' c E) > n c M) [f] Similarly, one can also ex- 



pect that 



> v. 



(A). 



in the flat A model, the assumed value of fl n 



6 An alternative explanation of this inequality follows from the results of sub- 
subsection 13.2.21 In the EdS model there is a large gravitational blueshift, so 
for the total redshift to be equal to that in the Milne model, a Doppler shift 
must be much larger. In other words, v e * must be much closer to the velocity 
of light than i> e M ' . 
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is smaller than unity; moreover, there is a late-time phase of accel- 
eration due to the cosmological constant. Interestingly, however, the 
value of De A) is greater than vi M> and much closer to ?4 B) . Appar- 
ently, for such a high redshift (and the adopted value f2 m = 0.24), 
the effect of initial deceleration of the Universe on recession veloci- 
ties strongly prevails over the opposing effect of its later acceleration. 
In fact this is not so surprising, since for f2 m = 0.24 the flat model 
starts to accelerate only for the redshift as small as about 0.5. 

On the other hand, for an arbitrarily high redshift there is always 
a range of sufficiently small values of f2 m , for which Uo A) < v& . 
We have checked that for f2 m = 0.24, the critical redshift for which 
«c A> = Uc M> is about 1.6. 



6 DISCUSSION 

Chodorowski (2007a) noted that in conformally flat coordinates in 
open, zero-A models, coordinate velocities of all galaxies are con- 
stant and subluminal. Lewis et al. (2007) pointed out that even in 
these coordinates, the proper velocitie of sufficiently distant galax- 
ies are superluminal. Lewis's et al. (2007) criticism was fully justi- 
fied, because coordinate velocities are meaningless. However, in this 
Paper we have shown that the proper velocity is not a real velocity of 
a galaxy, so the fact that the former can be superluminal does not im- 
ply the same for the latter. Indeed, basing on the definition of relative 
velocity given in this Paper we have demonstrated that velocities of 
all particles of the cosmic substratum are in fact not superluminal. 

In a recent paper, Faraoni (2010) studied the generation of the 
cosmological redshift by an instantaneous expansion of otherwise 
Minkowski spacetime. A source and an observer were assumed to 
be at rest relative to each other both at the time of emission and ob- 
servation of photons. It was also assumed that the universe expands 
only at a single instant of time tx and t e < tx < t . Faraoni (2010) 
showed that the induced redshift is given by the canonical formula 
for all FLRW universes, 1 + z = a(t )/a(t c ). On the basis of this 
finding and the fact that the source and the observer were initially 
and finally at rest, he deduced that in the studied case the cosmolog- 
ical redshift is entirely gravitational. This is in agreement with our 
results, as will be explained below. Equations J20 M2U are valid for 
any function a(t), not necessarily determined by the FL equations. 
The assumptions made by Faraoni imply that a(t e ) = a(t Q ) — 
and from Equation J2U one can see that, both at the time of emis- 
sion and observation, the Dopplerian component of the cosmological 
redshift is zero. [Alternatively, from Equation j20\ it follows that the 
recession velocity vanishes.] However, Faraoni's (2010) conclusion 
that the cosmological redshift is in general gravitational is incorrect. 
The quantity a vanishes only for closed FLRW models and only for 
a single instant of time. More importantly, our Universe is certainly 
expanding now, and whatever the exact values of cosmological pa- 
rameters, its expansion has been continuous in the past. Therefore, 
the cosmological redshift must be partly Dopplerian. 

BH9 transported the four-velocity of a distant galaxy along the 
null geodesic connecting the source and the observer. [The original 
idea dates back to Synge (1960) and Narlikar (1994)]. Instead, as the 
path of parallel transport we have chosen the geodesic connecting the 
source and the observer on the hypersurface of constant cosmic time. 
Both definitions are equally justified mathematically: both are intrin- 
sic (i.e. they do not depend on any coordinate system) and reduce to 

7 In this case, the proper distance is measured on the hypersurface of constant 
conformal time. 



the SR definition of relative velocity if the two observers are at the 
same point. The reason that as a path of transport we have chosen a 
geodesic on the hypersurface of constant cosmic time was to sepa- 
rate the effects of the curvature of the Universe from the effects of 
its evolution. We wanted to define the recession velocity at a well- 
specified instant of time. Transporting along a null geodesic provides 
vbh, which is a sort of time-average (over the time interval between 
the times of emission and observation) of our, instantaneous reces- 
sion velocity. Indeed, it can be easily verified that for small redshifts, 
hence also small time intervals, kbh = (v e + v )/2. As a conse- 
quence of their definition of relative velocity, BH9 found that in any 
FLRW cosmology and for any value of the redshift the cosmological 
redshift is related to the transported velocity vbh by the relativistic 
Doppler formula|f] From this fact they deduced that the cosmological 
redshift is entirely Dopplerian. Our interpretation of the redshift is 
different: in a non-empty Universe there is gravitational field, induc- 
ing a gravitational shift. Therefore, with an exception of the empty 
model, the origin of the cosmological redshift must be partly gravi- 
tational. Using our definition of relative velocity, we have thoroughly 
proven this conjecture in the present Paper. Thus, the difference be- 
tween the finding of BH9 and ours is yet another hint that vbh is not 
the recession velocity of a galaxy at a given instant of cosmic time. 
Using Equations l |37t and J4U it can be shown that for small redshifts, 
z = (/3 C + /3o)/2 + [(p e + /? )/2] 2 /2. Comparing this equation with 
Equation J4U we see that the relation between z and the effective ve- 
locity (/3 e + /9o)/2 is indeed Dopplerian. But this effective velocity is 
just vbh. 

Different definitions of relative velocity lead to different inter- 
pretations of the cosmological redshift. However, the definition of 
BH9 and ours imply the same very important conclusion: recession 
velocities of all galaxies inside our particle horizon are subluminal. 
This fact is perhaps more fundamental than the specific value of re- 
cession velocity for a given redshift. 



7 SUMMARY 

In this Paper we have decomposed the cosmological redshift into a 
Doppler shift and a 'rest', which we have interpreted as a gravita- 
tional shift. In order to extract the kinematic component from the cos- 
mological redshift, we had to define properly the recession velocity 
of a distant galaxy. This velocity is a relative velocity of the galaxy 
and the observer, and in GR one cannot directly compare vectors at 
widely separated points. First, one has to transport parallely the four- 
velocity of the distant galaxy to the observer. Next, the transported 
four-velocity is transformed to the local inertial coordinates of the 
observer and then the recession speed is typically extracted from the 
(transformed) radial component of the four-velocity. In general, par- 
allel transport depends on the chosen path, but geodesies are special 
curves in GR and they seem to be a very natural choice for the prob- 
lem under consideration. As the path, Synge (1960), Narlikar (1994) 
and BH9 adopted the null geodesic of photons emitted towards the 
observer. For reasons explained in Section[6] we have chosen a differ- 
ent path: the geodesic connecting the source with the observer on the 
hypersurface of constant cosmic time. BH9 found that in any FLRW 
cosmology and for any value of the redshift, the cosmological red- 
shift is related to the transported velocity, as defined by them, vbh, 
by the relativistic Doppler formula. From this fact they deduced that 

8 Which means that, in particular, for small redshifts the formula of BH9 does 
not reduce to the non-relativistic limit of Bondi (1947) 
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the cosmological redshift is entirely Dopplerian. However, from the 
Principle of Equivalence (PoE) it follows that the cosmological red- 
shift must be partly gravitational. The explanation of this apparent 
paradox is that kbh is not a recession velocity at a specific instant of 
time, but rather a sort of an effective, time-averaged velocity. 

In Section[3]we have presented analytical relations between the 
cosmological redshift and its Doppler component for two particularly 
simple cosmological models. In an empty universe, the cosmological 
redshift is entirely kinematic. This is expected, because in the empty 
model there is no gravitational field. In an Einstein-de Sitter universe, 
at the time of emission of photons the Doppler shift is greater that the 
total shift, so a gravitational shift is in fact a blueshift. At the time 
of observation, the situation is different: for large enough redshifts 
the gravitational component is dominant; in particular, even the par- 
ticle horizon is receding with a speed smaller than that of light. This 
marked difference in the magnitude of kinematic components of the 
total redshift at the time of emission and at the time of observation is 
due to the fact that all matter-only models have expanded faster in the 
past. 

Using only the Principle of Equivalence, in Section[4]we have in- 
dependently decomposed the (non-relativistic) cosmological redshift 
into a Doppler shift and a 'rest', and have shown that in this regime 
the 'rest' is induced by the Newtonian gravitational potential. The 
obtained result is in agreement with a second-order expansion of the 
exact formula derived here in the framework of GR (Sec. [3}, as well 
as with the classical formula of Bondi (1947) and equivalent findings 
of Peacock (1999) and Gr0n & Elgar0y (2007). 

Last but not least, in Section [5] we have critically revised a 
widespread opinion that the recession velocities of (sufficiently) dis- 
tant galaxies are superluminal. We have demonstrated that the com- 
monly used proper velocity is not a recession velocity, but merely 
a radial component of the galaxy's velocity four-vector (when ex- 
pressed in proper coordinates). As such, it can indeed attain super- 
luminal values. On the contrary, the definition of the recession ve- 
locity, presented in this Paper, naturally implies that the velocity is 
subluminal for any value of the galaxy's redshift. The definition of 
BH9 implies exactly the same conclusion. Basing on our definition 
we have also argued that the expansion of the Universe is subluminal 
even on superhorizon scales. We plan to study this issue more com- 
prehensively in a future work. 



Faraoni, V. 2010, GReGr, 42, 851 

Fraknoi, A., Morrison, D., Wolff, S. 2004, Voyages Through the Universe, 

(Belmont, CA: Thomson Brooks/Cole), 3rd ed. 
Francis, M. J., Barnes, L. A., James, J. B., Lewis, G. F. 2007, PASA, 24, 95 
Gr0n, 0., Elgar0y, 0. 2007, AmJPh, 75, 151 

Harrison, E. 2000, Cosmology, (Cambridge: Cambridge University Press), 
2nd ed. 

Kaufmann, W. J., Freedman, R. 1999, Universe, (W. H. Freeman, New York), 
5th ed. 

Lewis, G. F, Francis, M. J., Barnes, L. A., James, J. B. 2007, MNRAS, 381, 
L50 

Lewis, G. F, Francis, M. J., Barnes, L. A., Kwan, J., James, J. B. 2008, MN- 
RAS, 388, 960 
Milne, E. A. 1933, Zeitschrift fur Astrophysik, 6, 1 

Milne, E. A., McCrea, W. H. 1934, QuartJMath, 5, 73; reprinted in 2000 by 

GRG, 32, 1949 
Narlikar, J. V. 1994, AmJPh, 62, 903 

Padmanabhan, T. 1993, Structure Formation in the Universe, (Cambridge, UK: 

Cambridge University Press), p. 54 
Peacock, J. A. 1999, Cosmological Physics, (Cambridge, UK: Cambridge 

University Press), p. 87 
Peacock, J. A. 2008, preprint <arXiv:0809.4573l 

Peebles, P. J. E. 1993, Principles of Physical Cosmology, (Princeton, NJ: 
Princeton University Press), p. 94 

Rindler, W. 1977, Essential Relativity, (New York, Heidelberg, Berlin: 
Springer- Verlag), revised 2nd ed. 

Schutz, B. F. 1985, A First Course in General Relativity, (Cambridge, UK: 
Cambridge University Press) 

Seeds, M. A. 2007, Foundations of Astronomy, (Belmont, CA: Thomson 
Brooks/Cole), 10th ed. 

Synge, J. L. 1960, Relativity: The General Theory, (Amsterdam: North- 
Holland) 

Whiting, A. B. 2004, The Observatory, 124, 174 
Will, C. M. 2006, Living Reviews in Relativity, 9 



ACKNOWLEDGEMENTS 

This work was partially supported by the Polish Ministry of Science 
and Higher Education under grant N N203 0253 33, allocated for the 
period 2007-2010. 



REFERENCES 

Abramowicz, M. A., Bajtlik, S., Lasota J.-R. Moudens, A. 2007, AcA, 57, 139 
Abramowicz, M. A., Bajtlik, S., Lasota J.-R. Moudens, A. 2009, AcA, 59, 131 
Barnes, L. A., Francis, M. J., James, J. B., Lewis, G. F. 2006, MNRAS, 373, 
382 

Bondi, H. 1947, MNRAS, 107, 410 
Bunn, E. F, Hogg, D. W. 2009. AmJPh, 77, 688 (BH9) 
Chodorowski, M. J. 2007a, Old New Concepts Phys., 4, 15 
Chodorowski, M. J. 2007b, MNRAS, 378, 239 
Chodorowski, M. J. 2008, preprint (arXiv:0812.3972l 
Cook, R. J., Burns, M. S. 2009, AmJPh, 77, 59 
Davis, T. M., Lineweaver, C. H. 2001, AIPC, 555, 348 
Davis, T. M., Lineweaver, C. H. 2004, PASA, 21, 97 



© 0000 RAS, MNRAS 000, 000-000 



